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Fundamental mechanisms in premixed 
flame propagation via vortex-flame 
interactions - numerical simulations 


By T. Mantel 
1. Motivation and objectives 

During the past few years, direct numerical simulations (DNS) have been exten- 
sively used to study turbulent reacting flows in order to obtain a better understand- 
ing of the interaction between a turbulent flow field and a flame front, mainly for 
modeling purposes. We can cite different studies with increasing degree of com- 
plexity. Three-dimensional DNS of a decaying isotropic turbulence with chemical 
reactions (no heat release and constant density) have been carried out by Picart 
et al. (1988), Rutland et al. (1990), El Tahry et al. (1991). More recently, heat 
release and variable density have been taken into account by Poinsot et al. (1991), 
who have analyzed the interaction between a two-dimensional vortex pair and a 
premixed laminar flame. Haworth & Poinsot (1992) and Trouve & Poinsot (1994) 
focused their attention on the effects of the Lewis number by studying the inter- 
action between (respectively) a two- and three-dimensional decaying homogeneous 
turbulence and a premixed flame front. We can also note the work of Poinsot et 
al. (1994), who performed a two-dimensional interaction between a turbulent flame 
front and a cold wall. 

Due to the limitations of the available computers, all the studies cited above 
have been performed using a one-step irreversible reaction (Reactants — ► Products) 
to describe the chemistry occurring in a flame. Consequently, some conclusions 
of these studies remain questionable due to complex chemistry effects. This last 
remark is supported by the work of Baum et al. (1992), who have studied the 
interaction between a two-dimensional decaying isotropic turbulence and a stoi- 
chiometric hydrogen-air premixed flame using the 9 species, 19 reactions scheme of 
Miller et al. (1982). Despite the Lewis number (based on H 2 ) of their simulation 
being definitely less than unity, the local flamelet speed S„ decreases with increas- 
ing tangential strain rate, which corresponds to the behavior of a flame having a 
Lewis number greater than unity (Clavin 1985, Law 1988). The Lewis number is 
defined as the ratio of the thermal diffusivity (in the fresh gases) a u and of the 
molecular diffusivity of the limiting species T>: Le = a u /T>. Moreover, they find a 
better correlation between the tangential strain rate and S„ rather than the local 
curvature and S n . 

This contradicts partially the previous work of Haworth & Poinsot (1992), who 
studied the same configuration using a one-step Arrhenius chemistry. In that study, 
they observe that the curvature seems to be the determinant parameter controlling 
the local flame structure. They notice a strong correlation between the curvature 
and the local flame velocity 5 n for all the Lewis numbers investigated (Le = 0.8, 
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1.0, and 1.2). For the tangential strain rate, they found a significant correlation 
(with S n ) for Le = 1.0. 

This short discussion points out the need to verify the limit of validity of sim- 
ple mechanisms able to be implemented in DNS codes. Because reduced chemical 
mechanisms are available and represent a good compromise between (too) simple 
mechanisms and (unreachable) full chemical schemes, it is of primary interest to 
investigate the behaviors of some of these models on a configuration which exhibits 
the fundamental mechanisms occurring in premixed turbulent flames and for which 
experimental results are available. 

The interaction between vortex pairs or ring vortices and a premixed planar flame 
possesses the major features encountered in turbulent premixed combustion such as 
unsteadiness, stretch, curvature, Lewis number, radiative heat losses, and complex 
chemistry effects. This simplistic configuration has been previously studied both 
experimentally and numerically for different purposes by several authors (Poinsot 
et al. 1991, Rutland &; Ferziger 1991, Lee et ai 1993, Roberts et al 1992-1993, 
Samaniego et al. 1994-b, Driscoll et al. 1994). Some of these works have focused 
their attention on the effects of radiative heat losses and of the strain on the quench- 
ing of a premixed flame. The motivation of these studies was to propose an update 
of the premixed turbulent combustion diagrams to clarify the limit between the 
flamelet and non-flamelet combustion regimes (Poinsot et al. 1991, Roberts et al. 
1993). Other studies have extracted different statistics of curvature and orientation 
factors of turbulent premixed flame fronts (Lee et ai 1993). 

The goal of the present study is to assess numerically the ability of single-step 
and two-step chemical models to describe the main features encountered during the 
interaction between a two-dimensional vortex pair and a premixed laminar flame. 

This paper represents the second part of a joint experimental and numerical 
project concerning vortex flame interactions performed at the Center for Turbu- 
lence Research. This first part investigated by Samaniego et ai (1994-b) concerns 
the experimental aspect of this project. Thus, the configuration retained in our 
study corresponds to the experimental one of Samaniego et ai (1994-b). Briefly, it 
concerns the interaction between a two-dimensional vortex pair generated by acous- 
tic excitation and a V-shaped flame stabilized on a heated wire. In the experiment, 
imaging of the light emitted by the flame and smoke visualization of the flow field 
have been carried out to provide initial conditions to the simulation and eventually 
data to perform quantitative comparisons between the experiment and the simula- 
tion. Light emission imaging allows determination of the time history of the flame 
surface area, of the global heat release, and of the distribution of the heat release 
along the flame front. The characteristics of the vortex pair (circulation, position of 
the vortices with respect to the flame, distance between the center of the vortices) 
are obtained from smoke visualization. 

In the two-step mechanism, the reaction kinetics are represented by a first chain 
branching reaction A + X — ► 2X and a second chain termination reaction X -f X — + P 
(Zel’dovich 1948). This mechanism has been successfully used to analyze different 
features of premixed laminar flames (Linan 1971, Hocks et ai 1981, Seshadry and 
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Peters 1983). In particular, Seshadry and Peters (1983) investigated the response of 
a premixed laminar flame to an external strain. They pointed out the relative role 
of the Lewis number for the reactant A and the intermediate species X. In the case 
of a positive stretch, they show that the high diffusivity of the intermediate species 
contributes to a decrease of the reaction rate independently of the Lewis number of 
the reactant species. 

This paper presents the fundamental mechanisms occurring during vortex-flame 
interactions and the relative impact of the major parameters encountered in turbu- 
lent premixed flames and suspected of playing a role in quenching mechanism: 

- Influence of stretch is investigated by analyzing the contribution of curvature 
and tangential strain on the local structure of the flame. The effect of Lewis 
number on the flame response to a strained field is analyzed. 

- Radiative heat losses which are suspected to be partially or totally responsible 
for quenching (Poinsot et al. 1991, Roberts et al. 1993) are also investigated. 

- The effect of the diffusion of radicals is studied using a two-step mechanism 
in which an intermediate species is present. The parameters of the two-step 
mechanism are entirely determined from physical arguments. 

- Precise quantitative comparisons between the DNS and the experimental results 
of Samaniego et al. (1994-b) are performed. These comparisons concern the 
evolution of the minimum heat release rate found along the flame front during 
the interaction and the distribution of the heat release rate along the flame front. 


2. Accomplishments 

2.1 The mathematical model 


2.1.1 The conservation equations 

The DNS code has been developed at the Center for Turbulence Research fol- 
lowing the methodology of Lele (1992) and Poinsot & Lele (1992). The code fully 
resolves the compressible Navier-Stokes equations using a sixth order spatial scheme 
(Lele 1992) and a third order temporal scheme (Wray 1990). Due to the heat release 
and the resulting gas expansion, the boundary conditions of the computational do- 
main are inflow/outflow (Poinsot & Lele 1992). The transport equations solved in 
the DNS code can be written 


dp 

dt 
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where 
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(2.5) 


The transport equation for the species a is classically written: 


dpY Q 

dt 



( pUjYa ) = - 


djf 

dxj 


+ W a 


( 2 . 6 ) 


where is the diffusion of the species a flux modeled using a Fickian approxima- 
tion and w a the sink term due to chemical reaction. 

Here pE represents the toted energy, Q is the heat released by chemical reaction 
per unit mass of fresh gas, the subscript u corresponds to the unburnt gases. (Q 
ernd w a will be defined in detail later in this section for the one-step and two-step 
chemical models.) In the energy equation, the heat losses axe represented by a term 
linear in temperature and a heat loss parameter h which will be described later. 

The transport properties of the fluid are temperature dependent following a power 
law: 



where p represents the dynamic viscosity and b is a constant (here b = 0.76). 
The thermal conductivity A and the molecular diffusivities V a for the species are 
determined by assuming constant Prandtl and Schmidt numbers. 

2.1.2 The one-step chemical model 

In this model, the chemistry is described by a single step irreversible reaction: 


A(reactants) —> P(products) (2.8) 

The reaction rate of this reaction is expressed using a classical Arrhenius law. Fol- 
lowing the notations of Williams (1985) the reaction rate for the deficient species 
(subscript A) can be written 


ii>A = A pY a exp 



-e) \ 

1 — Qr(l — 0)/ 


(2.9) 


where 


0 —(T — T u )/(Tb — T u ) is the reduced temperature; 
a =(Ti> — T u )/Tft represents the heat release parameter; 

j3 =aE a / R°Ti, is the Zel’dovich number; 

A =B exp(—f3/a) is the pre-exponential factor. 


Fundamental mechanisms in premixed flame 


49 


Tb and E a represent respectively the temperature in the burnt gases and the acti- 
vation energy of the reaction, R° being the universal gas constant. 

The source term Q appearing in Eq. 2.3 is 

Q = (-A° f )w A (2.10) 

where (— A^) is the heat of reaction per unit mass of reactant. 

2.1.3 The two-step chemical model 

The two-step mechanism initially proposed by ZePdovich (1948) consists of a first 
order chain branching reaction and a second-order termination reaction: 

A + X — > 2X ( 2 . 11 ) 

X + X-^P (2.12) 

In the first reaction, two radicals X are created while one radical is consumed 

during the transformation of the reactant A into X. This initiation step is essential 
because it provides radicals to initiate chain reactions. This type of reaction is 
considered as thermo-neutral and has a high activation energy. Then, two radicals 
recombine to form the product P during the termination step (or chain breaking). 
This recombination step is highly exothermic and all the energy is release during 
this step. These approximations are consistent with the simplistic description that 
the more exothermic a reaction is, the smaller the activation energy. 

Thus, using this description and considering equal molecular weight 
(W = Wa = Wx), the reaction rates for the two reactions axe 

% ) (2.13) 

RR2 = ^P 2 Y£ (2.14) 


and for the two species, 


w A = -\ iP 2 y a y x ex p(-T3^T3^)) ( 2 - 15 ) 

wx = Ai p 2 Y A Yx ex p( - Y~^i^T 0 )) ~ 2 Az P 2y x (2.16) 

where 



In the energy equation, the source term due to chemical reactions is 

R 

Q = £(-A H° k )RR k 

k=l 


(2.17) 
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where A H® is the heat released by the k tk reaction. In the present case, we consider 
that the first reaction is thermo-neutral and all the heat is liberated during the 
second reaction. Thus Q becomes 


Q = (-A H° 2 )^p 2 Y> . 


(2.18) 


This mechanism has been utilized in various studies concerning premixed laminar 
flames. Linan (1971) proposed to represent a premixed laminar flame as a succession 
of layers in which different transformations occur. On the unburnt gas side, a 
thick preheat zone (relative to the laminar flame itself) is present. In this layer, 
only diffusion and convection phenomena take place and chemical reactions do not 
exist. Then, the chemistry of the hydrocarbons takes place in a thin layer in which 
radicals axe produced. These radicals are transported (by diffusion and convection 
towards the burnt gases and by diffusion towards the fresh gases) and recombine into 
products in a broader layer including the hydrocarbon consumption layer. Different 
regimes were encountered by Linan depending on the ratio of the frequency factors 
of the two reactions A 2 /A 1 . 

Hocks et al. (1981) investigated the quenching processes related to the interaction 
between a premixed laminar flame and a cold wall. The authors conclude that the 
two-step mechanism is able to describe the mechanism leading to flame quenching. 
It appears that the behavior of the flame close to the wall is strongly dependent 
on the maximum concentration of the intermediate species, which is essentially 
determined by the ratio A 2 /A 1 . 

Later, Seshadry & Peters (1983) studied the structure of a planar premixed lam- 
inar flame submitted to stretch. Considering a high activation energy for the first 
reaction, the authors derived an asymptotic expansion for the temperature. They 
found that the first order temperature can be expressed as a function of the stretch 
and the Lewis numbers for the reactant and the intermediate species: 



Le A ~ 1 
Le A 


+ (-a h* 2 ) 


1 ~ Le x 
Le x 



(2.19) 


The subscript 0 refers to the axial coordinate where Y\ is maximum, I is a func- 
tion always positive, K* is a non-dimensionalized stretch and (— AiT^ ) is the non- 
dimensionalized heat of reaction of the recombination step. The relation (2.19) 
points out the respective roles of the diffusivities of the reactant and of the interme- 
diate species. Considering only the first term on the RHS of Eq. 2.19, for positive 
stretch the temperature decreases (increases) for Le A > 1 (Le A < 1). For Le A — 1, 
the temperature remains constant equal to the zero order temperature, regardless 
of the value of the stretch. This recovers the classical conclusions of the role played 
by the Lewis number of the reactant on stretched flames (Clavin 1985, Law 1988). 

The second term on the RHS of (2.19) enhances the effects of diffusivity of the 
intermediate species on the dynamics of stretched flames. Since radicals are mostly 
very light species, they have high diffusivities leading to Lewis numbers significantly 
less than unity. Thus, in case of positive (negative) stretch, the diffusivity of the 
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intermediate species tends to decrease (increase) the temperature and then the local 
laminar flame speed. This result points out that even for Le^ = 1, the flame can 
be sensitive to the effect of stretch and exhibits locally a variation of the laminar 
flame speed. 

This brief description shows the ability of the two-step mechanism to describe 
fundamental phenomena occurring in premixed flames in different configurations. 
However, it also points out the appearance of new parameters such as the ratio of 
the frequency factors A 2 /A 1 , the Lewis number for the intermediate species, the 
activation energy and the heat released by each of the two-step mechanism. This 
represents an infinity of combinations between these parameters. The challenge, 
then, is to provide this model with a realistic set of parameters representative of 
the fundamental mechanisms encountered in the kinetics of hydrocarbons. 

To do so, we will focus our attention on the hydrogen-oxygen submechanism 
which is hierarchically the first submechanism (followed by the oxidation of carbon 
monoxyde) occurring in the chemistry of hydrocarbons (Glassman 1987, Westbrook 
and Dryer 1984). In this submechanism, the H atoms play a determining role 
because they provide a major source of radicals in the branching reactions of the 
oxidation of H 2 . Furthermore, the concentration of H radical directly affects the 
overall heat release and, consequently, the reaction rate (Westbrook and Dryer 
1984). In the H 2 — 0 2 submechanism, one of the most important chain branchings 
is 

H + 0 2 -> OH + 0 . (2.20) 

In the chain branching reaction (2.20), one H atom is consumed whereas radicals 
O and hydroxyl OH are produced and contribute to further branching reactions in 
which H radicals are created such as 

O + H 2 — ► OH + H 
OH + H 2 ^H 2 0 + H 
OH + OH -+H 2 O + O 

The reaction 2.20 is endothermic of 17 kcal/mol and has an activation energy of 
14.4 kcal/mol (Yu et ai 1994). 

Then, radicals are transformed to form product and liberate energy in chain 
breaking reaction. In the high temperature regime, a principal termination reaction 
is 

H + OH + M -> H 2 0 + M , (2.21) 

where M is a third body. During this step, radicals H and OH recombine to 
form products of combustion (here water). This reaction is highly exothermic (120 
kcal/mol) and has a zero activation energy. 

The rate coefficients for these two reactions are reported in Table 1. 

From these considerations, all the parameters appearing in the two-step mecha- 
nism 2.11 and 2.12 can be estimated from the rate coefficients reported in Table 1. 
The ratio of the frequency factors A 2 /Ai of reactions 2.11 and 2.12 can be deter- 
mined from the values of B\ and 1? 2 and from the temperature of the burnt gases. 
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k=BT n e~ 

E A /R° 

T 



B 

n 

E a 


reaction 

(cm 3 rnol“ 1 s” 1 ) 


(kJ/mol) 

ref 

ho 2 -oh+o 

8.3 10 13 

0 

60.3 

Yu et al. 1994 

h+oh+m-h 2 o+m 

1.6 10 22 

-2 

0 

Miller and Bowman 1984 


Table L Reacting mechanism, rate coefficients 


As we mentioned previously, the II atoms play a determinant role in the chemistry 
of hydrocarbons. Thus, we will identify the intermediate species X of the two-step 
mechanism to the H atom which leads to a Lewis number for X, Lex = 0.15. The 
activation energies of the two-step mechanism are those of the reactions 2.10 and 
2.11 reported in Table 1. 

2.1.4 Initialization of the two-step mechanism 

To initialize the one-dimensional laminar flame using the two-step mechanism, the 
assumption of quasisteady state for the intermediate species is considered. Thus 
from Eq. 2.16, the concentration of the radical is directly related to the concentra- 
tion of the species A and to the temperature: 




ft(i-e \ 

1 — a(l — 0)/ 


(2.22) 


We note that using Eq, 2.22 in Eq. 2.15 makes the two-step mechanism (under 
the assumption of quasi steady state for the intermediate species) reduce to the 
one-step second-order mechanism, 


A + A — ► P . 


(2.23) 


2.1.5 The radiative heat loss model 

The radiative heat losses are taken into account in the energy equation by a 
linear term in temperature. The heat loss parameter h appearing in Eq. 2.3 comes 
from the asymptotic analysis of Williams (1985-pp 271-276). He considers the 
asymptotic structure of a one-dimensional premixed laminar flame submitted to 
radiation or conduction to a wall. Considering constant properties for the fluid, 
Williams proposes 

(2 - 24) 

where / is a constant and Pr a Prandtl number, v u represents the molecular viscosity 
in the fresh gases and the adiabatic laminar flame speed. 
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2.1.6 Configuration 

The configuration retained in this study concerns the interaction between a two- 
dimensional vortex pair and a planar premixed laminar flame. This geometry cor- 
responds to the experimental device of Samaniego et al. (1994). In the experiment, 
a mixture of methane or propane and air is introduced in a vertical channel of 
square cross section. The flame is stabilized on a heated wire and exhibits a two- 
dimensional V-shape. On the left side of the duct, a vortex pair is generated by 
acoustic excitation and interacts with the flame (see Fig. 1). Then, during the main 
duration of the interaction, the flow keeps its two- dimensionality. The DNS of all the 
experimental domain cannot be considered because of the high computational time 
and memory requirement. Nevertheless, the second flame (on the right of Fig. 1) 
does not play any significant role during the early stages of the interaction and 
allows us to consider only the left side of the domain where the interaction occurs 
(Samaniego et al. 1994). Moreover, one can demonstrate that due to the Gallilean 
invariance of the Navier-Stokes equations, the problem is equivalent to considering 
a sub-domain related to the frame of reference of the flame and convected by the 
mean flow field. 

The characteristics of the vortices (circulation, size, distance with respect to the 
flame, angle of impingement) are obtained from the experiment and are reported 
in Table 2. Here, the Damkohler is defined by Da = s /V dql u and the laminar 
flame thickness by 6 f = a u j S\. Three vortex pairs having different circulations and 
sizes are studied. For each interaction, a methane-air flame and a propane-air flame 
axe investigated in order to study the effects of the Lewis number. The methane- 
air flame with an equivalence ratio of 0.55 exhibits a Lewis number equal to unity 
whereas the propane-air flame with an equivalence ratio of 0.50 has a Lewis number 
of 1.8. The parameters of the different cases analyzed in this study are reported in 
Table 2. 

Since the structure of the vortex has not been precisely determined in the exper- 
iment, it is difficult to locate our study in the diagram of flame- vortex interactions 
proposed by Poinsot et al (1991). Nevertheless, from the PIV measurements of 
Driscoll et al (1994), the maximum tangential velocity of the vortices u™ ax and 
the radius of the vortex core a can be estimated. Driscoll et al. determined that 

2 < U 0 1 ( 1 X /Vd < 6 and 1/2 < a/s < 1/6, where Vd is the displacement velocity of 
the vortex pair and s the distance separating the center of the vortices. Thus, we 
can roughly locate our study in the quenching region of the diagram represented in 
Fig. 2. 

The temperature jump across the flame front is set exactly equal to 5, leading to 
the heat release parameter a = 0.8. 

Vortex pairs can be easily generated numerically, and several authors have pro- 
posed analytical solutions that satisfy the Navier-Stokes equations. Among these 
solutions, one of the most used in numerical simulations are the Oseen vortex (Os- 
een 1911) and the vortex “hat” used by Rutland (1989). The expression for the 
circulation T, the vorticity u ) z , and the tangential velocity uq are expressed in Table 

3 in function of the radial coordinate r, the vortex strength and the core radius 




1 10 100 1000 10000 


FIGURE 2. Diagram of vortex-flame interactions 
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vortex 1 

vortex 2 
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case 

Ti/slh 

Re i 

r 2 /s£*/ 

Re 2 

s/Sf 

D/6, 

Vd/S l 

Da 


1 

3286 

4600 

3569 

5000 

36.5 

104.8 

14.9 

m 



5657 

8000 

6686 

9333 

25.7 

104.8 

38.1 

0.53 



10286 

14667 

11314 



104.8 

66.8 

0.3 


1 

3286 

4600 

3571 


40.5 

104.8 

13.5 

2.91 

1.8 

2 

5714 

8000 

6667 

9333 

28.6 

104.8 

34.5 

0.81 


3 

10476 

14667 

11428 

16000 

28.6 

104.8 

61.0 

0.45 


Table 2. Dimensionless parameters for each configuration 


vortex “hat” 

Oseen vortex 

r 2 r 2 

~~r 2 

T = 27r'I' exp( 

/tZ 

r = ^(l — exp( 2)) 

* r 2 r 2 

u>z ~ ~^2 ( 2 — ^2 ) ex P ( — 9^2 ) 

ijl T 7 * 2 

ij) r 2 

W: = 0 ex P( 0) 

7T<T^ <7^ 

2 


Table 3. Circulation T, vorticity u? 2 , and tangential velocity ue for the vortex “hat” 
and the Oseen vortex 

a. The distribution of these quantities for the Oseen vortex and for the vortex 
“hat” are schematically represented in Fig. 3. 

Because its compact structure allows us to avoid numerical difficulties related to 
the boundaries of the domain (and then to reduce the size of the computational 
box), the vortex “hat” has been utilized in several studies (Rutland fa Ferziger 
1991, Poinsot et ai 1991). However, the vortex “hat” exhibits a rapid decrease of 
the tangential velocity and an inversion of the sign of the vorticity for r = \/2a , 
which is not representative of real vortices. During the interaction with the flame, 
this opposite vorticity generates undesirable stretch, leading to an artificial increase 
of the flame length. To avoid these problems and to perform future quantitative 
comparison with the experimental results of Samaniego et ai (1994), the Oseen 
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radius 



radius 



FIGURE 3. Characteristics of the vortices a) circulation T, b) vorticity w : , c) 
tangential velocity ug\ : Oseen vortex; : “hat vortex”. 

vortex is utilized in all our simulations. 

For these two-dimensional simulations, the computational domain is typically 
composed by 300000 grid points (750 x 400). 

2.1.7 Remarks 

In the simulations, some parameters have to be treated carefully in order to 
perform quantitative comparisons with the experiment. The interaction between the 
vortex pair and the flame can be entirely defined by the following set of parameters: 

- the ratio between the displacement velocity of the vortex pair and the laminar 
flame speed Vd/Sl\ 

- the ratio between the distance separating the center of the vortices and the 
laminar flame thickness s/8f, 

- the distance D separating the vortex pair and the flame; 

- the temperature jump across the flame front represented by the heat release 
parameter a; 

- the activation energy of the chemical reaction(s) represented by the Zel’dovich 
number (s) 

- the ratio between the frequency factors A 2 /A 1 for the two-step mechanism. 
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Although it sets the laminar flame velocity, the frequency factor A for the one-step 
chemical model (and Ai for the two-step mechanism) is not a determining parameter 
in this study. Indeed, the laminar flame velocity can be chosen arbitrarily as long as 
the ratio Vd/Sl is appropriate and the effects related to the compressibility of the 
flow are avoided. The only restriction is to limit the Mach number, based on the 
maximum velocity present in the flow field, to be below a critical value (i.e. 0.2). 
This restriction has an important practical impact, because it considerably reduces 
the computational time of the simulations. 

As we have already mentioned, the experiment allows us to determine the individ- 
ual circulation of the vortices and the distance separating the center of the vortices, 
the radius of the viscous core a and the exact distribution of the tangential velocity 
being unknown. However, the information provided by the experiment is enough 
to give a good representation of the interaction. Vortex pairs propagate by mutual 
induction, and, for point vortices, the displacement velocity of the ensemble is de- 
fined by the circulation of the vortices T and the distance separating the vortices s 


(Prandtl &; Tietjens 1934): 



(2.25) 


Although in the experiment the vortex cores definitely have a finite radius, some 
simulations performed on simple cases show that as long as the circulation of each 
of the vortices is conserved and the radius of the viscous cores is restricted to 
0 < a/s < 1/2, the dynamics of the interaction and the strain field applied to the 
flame are the same. Thus, in all the simulations presented here, the ratio a/s is set 
to 1/3. 

In other words, for given circulation T and distance s, vortex flame interactions 
can be represented either by point vortices or by finite viscous core radius and lead 
to the same conclusions. This last point has a very important impact on the elabo- 
ration of combustion diagrams for vortex flame interactions proposed by Poinsot et 
al (1991) usually defined by u^ ax / Si and a/6j . Thus, using the characteristics of 
the Oseen vortex reported on Table 3, a vortex flame interaction can be represented 
in the diagram by an infinity of points related by the relation: 


n m *Hpr(±Y l 

Sl 2tt \S,J 


(2.26) 


where Re v = T/ z/ is the Reynolds number of the vortex. Here, we propose to use 
the parameters Vq /Sl and s/8f to classify vortex flame interactions in the diagram 
of Poinsot et al (1991). Note that, since the ratios u™ ax /Vo and a/s in Poinsot et 
al (1991) and Roberts et al (1993) axe not fax from unity, the conclusions of their 
study axe still valid and are not contradicted by our comments. 

2.1.8 Post processing 

The experimental results obtained by Samaniego et al (1994) concern essentially 
the CH emission of the flame during the interaction. The CH emission of the flame 
is directly related to the heat release rate of the flame, as shown by the numerical 
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FIGURE 4. Evolution in time of the reduced flame length £* ; Case 1: : 

DNS, o : exp ; Case 2: : DNS, + : exp ; Case 3: : DNS, □ : exp 

analysis of Samaniego (1994). Samaniego et al relate the light emission of a 
strained laminar methane- air flame Ich 4 and the light emission of a strained laminar 
propane-air flame Iq 3 h 8 to the heat release rate RR according to the following power 
laws: 


RR = leu] (2.27) 

RR = &s (2.28) 

Thus, quantitative comparisons will be performed accurately both for methane-air 
and propane-air flames. 

In this section, local characteristics of the flame such as the tangential strain, the 
curvature, the normal towards the fresh gases, and the flame length E are obtained 
at the flame front location Ya/Ya u = 0.2 where Ya represents the mass fraction of 
the reactant A. 

2.2 Results 

In the first part of this section, the results are presented for simulations performed 
with the one- step chemical model. In this model, the heat release parameter a is 
set to 0.8 and the Zel’dovich number /? to 8.0, corresponding to an activation energy 
of 30 kcal/mol. 

Before discussing in more detail the physical effects affecting the local structure 
of the flame during the interaction, the evolution in time of the reduced flame length 
£* = £/£o for Le = 1.0 and the three interactions shown in Table 2 are presented 
in Fig. 4. 

The comparison with the experimental results of Samaniego et al. (1994) points 
out that the dynamics of the interactions is very well reproduced by the simula- 
tion. This enhances the reasonable approximations made concerning the internal 
structure of the vortices and their initialization. 
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FIGURE 5. Evolution in time of the minimum heat release rate, a) Le = 1.0 ; b) 

Le = 1.8 Case 1: : DNS, o : exp; Case 2: : DNS, + : exp; Case 3: 

: DNS, □ : exp 

Fig. 5 describes the time evolution of the minimum of the heat release rate along 
the flame front for the different interactions presented in Table 2. For Le = 1.0 
(Fig. 5-a), the agreement is poor between the simulations and the experiment even 
if the sensitivity to the Damkohler number is qualitatively represented by the DNS. 
For Le = 1.8 (corresponding to the propane-air flame), the simulations seem to well 
reproduce the decrease of the minimum heat release rate for the three interactions, 
especially those corresponding to low Damkohler numbers. 

The same comments can be made concerning the distribution of the heat release 
rate along the flame front where the simulations greatly overestimate the heat release 
for Le = 1.0 ( see Fig. 6), whereas the agreement is notably better for Le = 1.8 (see 
Fig - 7) - 

For Le = 1.8, the Lewis number effect alone can explain the decrease of the heat 
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FIGURE 6. Distribution of the heat release rate along the flame front. Le = 1.0 ; 
a): case 1, b): case 3 ; : DNS, : exp 

release rate during the interaction and seems to be the parameter controlling the 
local structure of the flame. However, the simulations presented here are performed 
for an adiabatic flame using a one-step chemical model. Thus, phenomena such 
as radiative heat losses and complex chemistry effects are not taken into account 
and could have a significant impact on the interaction, which could explain the 
discrepancies found between the DNS and the experiment, especially for Le = 1.0. 

In the following sections, the contributions of the major physical effects suspected 
of playing a role in the mechanism of extinction (stretch, heat losses, complex 
chemistry) and responsible for the decrease of the heat release rate observed during 
the interaction will be analyzed in detail. 

2.2.1 Analysis of the stretch 

In the configuration studied here, the flame is subjected to curvature and straining 
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FIGURE 7. Distribution of the heat release rate along the flame front 
Le = 1.8 ; a): case 1, b): case 3 ; : DNS, : exp 


effects, causing a modification of the flame length. These changes imposed on the 
flame can be estimated from the flame stretch «. From Williams (1985), k is defined 
by the rate of change of a Lagrangian flame surface element A: 


_ ^dA 

A dt 


(2.29) 


This expression can be rewritten in terms of flame stretch and flame curvature 
(Matalon 1983, Candel h Poinsot 1990): 

K = nn: Vw + V.w (2.30) 


where n is the unit vector normal to the flame surface pointing towards the fresh 
gases: 


n = — i 


Vc 
Vc | 


(2.31) 
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and w is the local velocity of the flame surface. In indicial notation, 
nn : Vw = niTijdw,/dxj. The local flame front velocity can be decomposed in 
a convective velocity in the fresh gases u and a displacement speed S d (see Poinsot 
et al. 1992). Thus, expression 2.30 becomes: 

rr S d 

« = V t .u — — (2.32) 

where V ( .u represents the strain contribution and S d /R the curvature term (7 Z is 
the local radius of curvature of the flame front). 

As it has been shown by Poinsot et al. (1992), the displacement speed S d is 
not equal to the laminar flame speed Si and can be significantly different. The 
displacement speed can be exactly determined by solving the equation for an iso- 
contour c = c 0 (which is nothing but the G-equation described in Kerstein et al. 
1988) as it has been proposed by Trouve & Poinsot (1994). The equation for an 
iso-scalar surface c = cq is written: 


dc 

dt 


+ w.Vc = 0 


(2.33) 


Utilizing Eq. 2.30 in Eq. 2.32 and noting that S d = w.n — u.n, an expression for 
S d can be derived: 


(2.34) 


Fig. 8 represents the relative contribution of strain and curvature on the stretch 
along the flame front for case 3 of Table 2 at the end of the interaction (t = 3.3ms). 

One sees that the maximum stretch does not appear on the trajectory described 
by the ensemble as it could be expected. The maximum of stretch occurs on both 
sides of this trajectory corresponding to each of the vortices. This feature was also 
observed by Driscoll et al. (1994) in their PIV measurement s . 

As expected due to the size of the vortex pair compared to the flame thickness, 
the contribution of the curvature on the global stretch is very weak (about a few 
percent) compared to the contribution of the tangential strain, especially in regions 
where a significant decrease of the heat release rate is observed. This corroborates 
the conclusions of Driscoll et al. (1994), who performed a measurement of the 
velocity field during vortex-flame interaction using a PIV technique. We can also 
see that the levels of strain rate reached during the interaction sire extremely high 
compared to the extinction strain rate measured for steady counterflowing premixed 
flames (Law et al. 1986, Tsuji & Yamaoka 1982) and for the vortex-flame interaction 
of Driscoll et al. (1994). In Tsuji & Yamaoka (1982), the extinction strain rate of 
a methane- air flame with an equivalence ratio of 0.52 is only 42s -1 . During the 
interaction of ring vortices and a premixed methane-air flame with an equivalence 
ratio of 0.55, the extinction strain rate measured by Driscoll et al. (1994) is 35s -1 
and corresponds to a Karlovitz number of 0.12. In the present interactions, the 
Karlovitz number reaches higher values up to 3.0 (see Fig. 9). 


II 
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Apparently, the configuration retained in our study exhibits a very resistant flame 
compared to the steady situations of Tsuji & Yamaoka (1982) and particularly 
compared to the ring vortices-flame interaction of Driscoll et al. (1994). In this last 
study, local quenching is clearly observed during the interaction. On the contrary in 
our study (both in the experiment and in the simulations), all the cases investigated 
here are far from the quenching limit encountered by Roberts et ai (1993) and 
by Driscoll et al. (1994). In all our simulations, the minimum heat release rate 
never decreases under 30% of the heat release rate corresponding to an unstrained 
laminar flame (see Fig. 5). Moreover, we observe that it is more difficult to quench a 
methane-air flame having a Lewis number of unity rather than a propane-air flame 
with Le = 1.8. These conclusions axe in accordance with the classical analytical 
results concerning the theory of stretched flames (Clavin Williams 1982, Clavin 
1985, Law 1988) and also with the experimental results of Law et al. (1986), but 
are in contradiction with the observations of Roberts et al. (1993). 

For Le = 1.0, we observe a good correlation between the strain rate and the 
local flame speed with a slight negative slope 1 (see Fig. 10-a). This is due to the 
compression of the reaction zone by the local extensive strain. This observation 
is consistent with the results of Haworth & Poinsot (1992) and is illustrated in 
Fig. 11, which shows the scatter plot of the thermal flame thickness 6^ = St/S^ 13 
(based on the maximum temperature gradient) versus the local flame speed. A very 
good correlation is found between 6^ and S£ and clearly shows a compression of 
the flame (down to 70% for the higher strain rates) corresponding to low values of 
S£. Conversely, a thickening of the flame (up to 10%) corresponds to higher flame 
velocities and regions of compression of the flame. 

For Le = 1.8 (see Fig.lO-b), the correlation between S£ and Ka is less obvious 
even if we observe a clear decrease of the local flamelet velocity for positive stretch 
as it can be expected for strained laminar flame having a Lewis number greater 
than unity. 

2.2.2 Effect of the radiative heat losses 

Radiative heat losses occurring in the burnt gases of a premixed laminar flame 
produce a natural decay of the temperature and have been suspected of being par- 
tially or totally responsible for quenching during vortex-flame interaction (Poinsot 
et al. 1991, Roberts et al. 1993). In their numerical study (and using a one-step 
mechanism for the chemistry), Poinsot et al. (1991) concluded that stretch alone 
cannot be responsible for quenching. Local extinction occurs only when stretch and 
radiative heat losses present in the burnt gases are combined. The authors observed 
that during the quenching process, pockets of fresh gases can be present in the burnt 
gases and cannot reignite due to a too low temperature. However, in order to ob- 
tain such behavior, very high (and unrealistic) heat losses have been imposed to the 
flame. In comparison, we will see that the heat losses measured by Roberts et al. 
(1993) and by Samaniego et al. (1994) are lower by an order of magnitude than in 

1 The local flame speed is calculated by integrating the reaction rate in a direction normal to the 

flame front: Sr = — b — f wdn 

u PuYa u j —oo 
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FIGURE 8. Contribution of tangential strain and curvature on stretch, Case 3 
le = 10 ; : total stretch, : curvature, : tangential strain 



time [ms] 


FIGURE 9. Evolution in time of the Karlovitz number for Le — 1.0 

: case 1, : case 2, : case 3 

Poinsot et ah (1991). Since in their simulations a single-step mechanism has been 
used, it is still not clear if quenching in real flames is due to combined stretch and 
heat loss effects, or to combined stretch and complex chemistry effects. 

Here, in order to quantify the impact of realistic heat losses on the distribution of 
the reaction rate along the flame front, a simulation corresponding to the interaction 
2 of Table 2 for Le = 1.0 is performed. Different levels of heat loss corresponding 
to the experiments of Roberts et ah (1993) and Samaniego et ah (1994) and to 
the numerical study of Poinsot et ah (1991) are investigated. Before discussing the 
results of these simulations, the quenching limit of a one-dimensional non-adiabatic 
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Karlovttz number 

FIGURE 10. Scatter plot of local flame speed versus Karlovitz number, 
a): Le = 1.0 ; b): Le = 1.8 


flame is investigated. This analysis should allow us to recover the asymptotic solu- 
tion of Williams (1985), who proposes an analytical relation between the laminar 
flame velocity and the heat loss parameter. This relation is expressed (Williams 
1985 p. 275): 

£ 2 = exp (“) (2.34) 

where £ = Sl/ S% d . From this relation, the quenching limit corresponds to / = 1/e 
leading to a reduced laminar burning velocity £ = 1 /y/e. 

Fig. 12 shows the evolution of the burning velocity ratio £ versus the heat loss pa- 
rameter l given by the simulation and by the relation 2.34 for two different Zel’dovich 
numbers (/? = 8 and j3 = 16). A very good agreement is found between the simu- 
lations and the asymptotic theory of Williams, especially for very large activation 
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local flame thickness 

FIGURE 11. Scatter plot of local flame speed versus the reduced local flame 
thickness 



FIGURE 12. Burning velocity ratio versus heat loss parameter ; asymptotic 

theory, o : /? = 8.0, + : /? = 16.0 

energy. 

In the simulation corresponding to the case 2 for Le = 1.0, various values for the 
heat loss parameter l are investigated. The case l — 0 corresponds to the adiabatic 
case, / = 0.008 and / = 0.027 to the heat losses estimated from the experiment 
of Samaniego et ai (1994) and from Roberts et ai (1993). The case / = 0.3 
represents the heat losses used by Poinsot et al (1991). According to Fig. 12, only 
a negligible effect of radiative heat losses is expected to be found in Samaniego et 
al (1994) and in Roberts et al (1992). Nevertheless, the effect should be much 
more pronounced for the heat losses of Poinsot et al (1991) which are close to the 
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FIGURE 13. Evolution in time of the minimum heat release rate, case 2 Le — 1.0 
; adiabatic, : l = 0.008, : l = 0.027, : / = 0.3. 

quenching limit observed for a one- dimensional unstrained non- adiabatic laminar 
flame. Fig. 13 represents the time history of the minimum heat release encountered 
along the flame front. As expected, the heat losses found in the two experiments 
have a negligible effect on the evolution of the heat release rate. On the other hand, 
very intense heat losses (more representative of conductive heat losses to a wall) can 
lead to a dramatic decrease of the heat release and to probable local quenching even 
if the computation does not go sufficiently far in the interaction to clearly show it. 

These simulations are performed using a one-step chemical model and do not 
take into account combined effects of radiative heat losses and complex chemistry. 
These combined effects can have a significant influence on the flame structure when 
the flame is close to its flammability limits, as illustrated by the numerical study 
of Egolfopoulos (1994). The author investigates the response of a lean methane- air 
flame submitted to stretch and radiative heat losses using the GRI (Gaz Research 
Institute) mechanism. For a laminar methane-air flame with an equivalence ratio 
of 0.55, Fig. 14 shows the evolution of the reduced heat release QlQ% d as a function 
of the tangential strain. Results for the adiabatic and non-adiabatic cases are 
presented. We can observe that the heat release rate is strongly influenced by the 
heat losses for the unstrained laminar flame and decreases down to 0.83. For high 
values of the strain (typically those encountered in our configuration), the effect 
of the radiative heat losses compared to the effect of the strain on the decrease of 
QIQ$ d is much weaker. 

2,2.3 Role of an intermediate species in strained laminar flames 

The chemistry of simple hydrocarbons such as methane can be described by full 
mechanisms available in the literature (Westbrook & Dryer 1984, Miller Bowman 
1989). The coupling of these mechanisms with a DNS code has been performed 
by Baum et al. (1992) but leads to extremely high computational time. As an 
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stretch [r l J 


FIGURE 14. variation of the heat release rate versus stretch for a one- dimensional 

premixed strained flame (from Egolfopoulos 1994) : adiabatic , : non- 

adiabatic 


alternate possibility, reduced chemical schemes for the combustion of methane-air 
flames (Peters & Williams 1987) and for propane-air flames as well (Kennel et al. 
1990) exist and can be used in DNS codes. 

Here, a two-step mechanism described by a first chain branching step A+X — ► 2X 
and a second chain terminating step X + X — ► P is utilized to analyze the structure 
of the flame during its interaction with a vortex pair. 

In the asymptotic analysis of a premixed laminar flame submitted to stretch, 
Seshadry & Peters 1983 have put into evidence the role played by the radical species 
on the flame structure (see §2.1.3). 

Since the Lewis number for the intermediate species is less than unity (due to 
the high diffusivity of the radicals), positive stretch will generate a high diffusion of 
radicals out of the reaction zone. Thus, the concentration of radicals in the reaction 
zone will decrease, leading to a lower heat release rate (because of the quadratic 
dependence on Yx of the heat release rate, see Eq. 2.14). This behavior is well 
reproduced by the asymptotic expression for the temperature given by Eq. 2.19. 

From the discussion presented in §2.1.3, the set of parameters for the two-step 
mechanism can entirely be determined by considering the data presented in Table 
1 and the temperature of the burnt gases (in the experiment T& = 1500A'). This 
leads to A 2 /A 1 = 0.53, = 4.83, /?2 = 0, and Lex = 0.15 (here using the Lewis 

number for the H atom). 

Fig. 15 shows the evolution in time of the minimum heat release rate during the 
interactions for Le = 1.0 and Le = 1.8 for the cases 1 and 3 of Table 2. Here, the 
decrease of Q is very well reproduced by the two-step mechanism, regardless of the 
Damkohler number. This is particularly interesting for the methane-air flame cases 
for which the one-step mechanism poorly describes the decrease of Q (see Fig. 5-a). 
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time [ms] 

FIGURE 15. Evolution in time of the minimum heat release rate, a) Le = 1.0 b) 

Le = 1.8 Case 1: : DNS, o : exp; Case 2: : DNS, + : exp; Case 3: 

: DNS, □ : exp. 

The coherent behavior of the two-step mechanism is confirmed by the distribu- 
tion of Q along the flame front (see Figs. 16 and 17). For Le = 1.0, the agreement 
between the DNS and the experiment is excellent and is greatly improved com- 
pared to the results obtained with the one-step mechanism (see Fig. 6). A minor 
quantitative discrepancy exists in regions of the flame where curvature effects are 
non- negligible although the tendencies axe well reproduced. 

The same comment can be made for Le — 1.8, where a good agreement is observed 
for the two Damkohler numbers studied here (Figs. 15-b and 17). However, the 
two-step mechanism does not significantly improve the results, which are already 
satisfactory using the one-step model (see Fig. 7). 

Fig. 18 shows the scatter plot of the local flame speed versus the Karlovitz number. 
Regions of the flame submitted to positive stretch exhibits a decrease of the local 
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curvilinear abscissa [mm] 

FIGURE 16. Distribution of the heat release rate along the flame front. Le — 1.0; 
a): case 1, b): case 3; : DNS, : exp. 

flame speed as it is observed for stretched flames having a Lewis number greater 
than unity. As a consequence, a flame having a Lewis number of unity described 
with a two-step mechanism behaves like a laminar flame having an apparent Lewis 
number greater than unity. This particular point provides evidence for the existence 
of a critical Lewis number Le c smaller than unity for which moderate stretch does 
not affect the flame structure. 

According to these observations, the concentration of radicals seems to be the key 
parameter which controls the local flame structure and quenching mechanism. This 
has been shown by the asymptotic analysis on the structure of strained premixed 
flames by Seshadry & Peters (1983) and is confirmed in the present study. Moreover, 
the quenching of a premixed laminar flame propagating toward a cold wall is also 
strongly dependent on the concentration of radicals (Hocks et ai 1981). 

However, it is difficult to separate which parameter of the two-step mechanism 
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FIGURE 17. Distribution of the heat release rate along the flame front. Le — 1.8; 
a): case 1, b): case 3; : DNS, : exp. 

is predominant (ratio A 2 /A 1 , Lewis number for the intermediate species, activation 
energy for the first and second step). For instance, to describe the flame quenching 
at a wall, Hocks et al (1981) have considered a low recombination regime, a high 
activation energy for the first reaction, and equal diffusivities of heat and mass 
(A 2 /A 1 = 5 1(T 4 , /?! = 8.5, #2 = 0, Lt x = 1.0). This choice has been made 
arbitrarily by the authors and justified by comparing the concentration of radicals 
given by the two-step mechanism to the concentration of intermediate species such 
as OH, H, and O given by complex chemistry calculations of the same configuration. 
Here, the choice of the set of parameters is based on physical arguments (justified 
in §2.1.3), which seems to be reasonable if we consider the agreement between the 
DNS and the experiment for all the cases studied here. The two-step mechanism 
shows its ability to describe vortex-flame interactions on a large range of Damkohler 
number (0.3 to 2.9). 
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Karlovitz Dumber 


FIGURE 18. Scatter plot of local flame speed versus the Karlovitz number. 

Conclusion 

Direct numerical simulations of flame vortex interactions are performed using 
one-step and two-step chemical models. 

Simulations performed with the one-step mechanism exhibit a strong disagree- 
ment between the DNS and the experiment for all the interactions with Le = 1.0. 
The simulations greatly overestimate the distribution of the heat release rate along 
the flame front during the interaction. On the contrary, for Le = 1.8, the agree- 
ment between the simulations and the experimental results is satisfactory for all the 
interactions corresponding to large range of Damkohler numbers (0.3 to 2.9). 

- The analysis of the stretch shows a negligible contribution of the curvature due 
to the large size of the vortex pair compared to the laminar flame thickness. The 
tangential strain generated by the vortices is responsible for the stretch felt by 
the flame. 

- Radiative heat losses representative of those encountered in the experiment have 
no effect on the local flame structure. However, tremendous heat losses (more 
representative of heat losses by conduction to a cold wall) have a dramatic effect 
on the heat release and lead to local quenching. 

- For the two-step mechanism, the set of parameters is entirely determined by con- 
sidering the rate coefficients of the chain branching reaction H + O 2 — ♦ OH -f O 
and the terminating reaction H+OH+M — ► H 2 O+M, which are the most impor- 
tant chain branching and chain breaking reactions of the H 2 — O 2 submech ani sm 
occurring in the chemistry of hydrocarbons. All the simulations performed with 
this model (for both Le = 1.0 and Le = 1.8) lead to very good agreement with 
the experimental results. During the interaction between the flame and the vor- 
tex pair, the concentration of radicals seems to be the key parameter controlling 
the local structure and the quenching mechanism. 
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As a consequence, it appears that the two-step mechanism with the set of param- 
eters chosen here is sufficient to describe all the interactions performed experimen- 
tally by Samaniego et al. (1994-b). 

Since most of the phenomena encountered in these interactions are also present 
in turbulent premixed flames, it is tempting to extrapolate these conclusions to a 
general description of turbulent premixed flames. 
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